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SECTION  I 
INTRODUCTION 


One  of  the  primary  considerations  in  the  simulation  of  an  isolated 
synchronous  generator  is  that  the  system  model  be  compatable  with 
loading  networks  of  many  different  varieties  and  yet  that  the  introduction 
of  changes  in  generator  parameters  or  loading  network  be  simply 
implemented.  A further  consideration  is  that  machine  parameters  be 
easily  identifiable  in  terms  of  the  results  of  physical  measurements. 

That  these  considerations  are  in  conflict  is  recognizable  by  the  fact  that 
current  standard  machine  practices  give  rise  to  measurements  that 
normally  identify  parameters  in  direct- quadrature  configuration,  while 
the  physical  machine  outputs  are  in  the  three  phase  configuration.  For 
many  loading  network  conditions,  a spatial  transformation  such  as 
Parks  or  Blondel  transformation  can  be  applied  to  the  loads  such  that 
loads  compatable  with  the  direct-quadrature  model  are  obtained.  How- 
ever, the  introduction  of  various  nonlinear  load  devices,  such  as 
rectifiers,  dictates  that  both  the  machine  representation  and  the  load 
n^'twork  representation  be  in  a normal  three  phase  formulation.  An 
additional  consideration  suggesting  compliance  with  three  phase 
representation  is  the  difficulty  in  inclusion  of  nonlinear  effects,  such 
as  saturation,  into  the  transformed  model. 

This  study  is  an  investigation  of  the  required  formulation  of 
a synchronous  machine  time  description  such  that  properties  previously 
discussed  are  incorporated  into  that  formulation  and  such  that  the 
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computation  be  time  efficient.  Also  addressed  in  this  study  is  the 
inclusion  of  an  effective  saturation  property  into  the  Sceptre  formulation 
such  that  a more  realistic,  a more  accurate  nonlinear  machine  model  is 
produced. 

SECTION  II 

LINEARIZED  MACHINE  REPRESENTATION 

The  electrical  schematic  of  a three  phase  machine  and  its  circuit 
equivalent  is  shown  in  figure  1.1.  The  equations  of  motion  for  this 
representation  are 

^ L ) + R i(t)  1.  1 

where  i(t)  = T i^(t),  i^(t),  i^(t),  i^lt),  iq(t) 

The  order  of  this  system  with  resistive  loading  is  considered 
to  be  of  sixth  order.  The  only  additional  increase  in  system  order 
necessary  could  be  through  inclusion  of  equivalent  field  driver 
circuitry,  the  possibility  of  which  is  included  into  the  resulting  simu- 
lation. The  i{j(t),  iq(t)  currents  represent  equivalent  direct  and 
quadrature  damper  currents  which  if  flowing  through  their  represen- 
tative circuits  would  produce  the  effects  of  the  machine  damper 
structure  or  amortisseur  windings  and  other  circulating  rotor  currents. 
Under  steady  state  operating  conditions,  the  damper  currents  are 
zero.  However,  with  the  introduction  of  a load  disturbance  the 
resultant  air  gap  flux  field  moves  relative  to  the  damper  windings, 
resulting  in  a damping  effect  on  the  electrical  transients  introduced 
into  the  machine  windings. 
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The  equivalent  resistance  in  each  circuit  is  represented  by 
the  six  by  six  diagonal  matrix, 

R = diag  R^,  Rj^^^  Rj^^  ^ 

The  v(t)  vector  is 

v(t)  =:  [va(t),  vt,(t),  Vc(t),  Vf(t),  0,  0 1 1.3 

The  phase  voltages  v^(t),  Vj^(t),  and  v^,(t)  are  voltages  which  result 
from  the  interaction  of  the  load  network  being  driven  and  the  generator 
Vf(t)  is  the  field  voltage  supply  and  for  the  simplest  model  is  assumed 
to  be  a source  with  high  internal  impedance  operating  at  constant 
current  conditions. 

L 18  a rotor  position  varient  inductive  matrix,  the  accepted 
standard  form  of  which  is  given  in  equation  1.4.  If  the  angular 
mechanical  frequency  of  the  machine  is  considered  to  be  a constant 
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(he  machine  equations  can  be  rewritten  in  the  forn'i 
L + L i(0  + R i(t)  = vtO 

whe  re 


a 


b 


L = 


-2u)La2  sin  2ail 
-2(i)La2  sin  (2oit-  2n/3) 
-2u)La2  sin  (2ujt  + 2n/3) 
-tJiMaf  sin  tut 
-a'Maf  sin  ii)t 
-xMaq  cos  ixt 


-2a' La 2 sin  (2x’t  - 2tt/3) 
-2aLa2  sin  (2ujt  f 2iT/3i 
-2irLa2  sin  2a’t 
-a'Maf  sin  (at  - Z.n/3) 
-a'Maf  sin  (at  - 2rr/3) 
-xMaq  cos  (at  - 2tt/31 


f 

-aMaf  sin  at 
-aMaf  sin  (at  - 2tt/3) 
-a'Maf  sin  (at  + 2n/3) 
0 
0 
0 


d 

-a'Maf  sin  at 
-aMaf  sin  (at  - 2tt/3) 
-a’Maf  sin  (at  + 2tt/3( 
0 
0 
0 


c 

-2aLa2  sin  (2a't  + Zn/3) 
-2aLa2  sin  2at 
-2aLa2  sin  (2a’t  - 2f/3) 
-aMaf  sin  (at  + 2r/3) 
-'oMaf  sin  (at  + 2n/3) 
-xMaq  cos  (at  + Z'r/i) 


q 

-OjMaq  cos  at 
-a'Maq  cos  (at  - 2tt/3) 
-a'Maq  cos  (at  + 2r/3t 
0 
0 
0 


This  assumption  of  constant  frequency  is  only  partially  valid  for 
most  airborn  auxiliary  power  units  - its  relaxation  in  the  resulting 
simulation  model  does  not  produce  traumatic  changes  in  formulation. 

If  the  prime  mover  model  is  available  the  torque  equation  must  be  used 
along  with  the  voltage  equation  because  the  rotor  position  becomes  a 
nonperiodic  function  of  time.  However,  in  general,  the  synchronous 
generator  has  very  little  variation  in  frequency  for  severe  load 
variations  and  (he  torque  equations  are  not  included  in  this  model. 


- S- 


The  resulting  system  of  equations  can  be  interpreted  by  examination 

of  the  first  phase  equation  which  takes  the  form 
r 

L^ao  La2  cos  Zuut 

r die  dif 

^ L^so  + ^a2  cos  (2(jjt  + 2tt/3)  J cos  tut  — j. 

did  . dig  o T ■ :>  ■ 

+ cos  o)t  + Mgq 

2(uL_^2  sin  (2u)t  - 2tt/3)  i^^  - 2(u  sin  (2tut  + 2ir/3)  i ^ 

- Mg^f  tu  sin  tut  i£  - tu  sin  tut  i^  - Mgq  tu  cos  ti)t  iq 

+ Ra  = '^a(t)  1.7 


aia  [-  -1  uiu 

dt  ^so  + l-'a2  ~ ^ dt 


The  equivalent  circuit  is  thus  seen  to  consist  of  self  inductance  with 
coupling  mutuals  plus  a sequence  of  dependent  sources  which  are 
functions  of  other  currents  inside  the  machine.  In  addition,  there 
exists  terms  including  the  winding  resistance  and  a time  variable 
resistance  in  the  phase  loop.  If  all  circuits  are  treated  in  the  fore- 
going dependent  voltage  source  manner  the  general  circuit  of  figure 
1.2  is  produced.  Furthermore,  the  form  of  most  of  the  dependent 
sources  is  such  that  they  are  representable  in  Sceptre  in  a simple 
equation  function,  A Sceptre  program  for  this  portion  of  the 
generator  model  is  listed  in  Appendix  A. 


SECTION  III 

MACHINE  PARAMETERS  IN  DIRECT  PHASE  FORMULATION 


In  general,  machine  parameters  are  obtained  experimentally 
in  quantities  directly  relatable  to  DQO  variables.  For  the  purpose  of 
completeness,  the  following  transformation  properties  from  DQO  to 
direct  parameters  is  presented. 

Consider  the  machine  equations  in  DQO  formulation  as  given 

in  [l]. 

v(t)  = R i(t)  + x(t)  + i((t),  x(t)  = L i(t)  1.8 

where 


n 

'"d 

■^D 

Q 

■‘Q 

^Q 

'"o 

^O 

^O 

i = 

^F 

X = 

Xf 

O 

'•KD 

^KD 

o 

^KQ  _ 

■^KQ  _ 

^DD 

O 

O 

^ DF 

^KDD 

O “ 

O 

a 

a 

O 

O 

O 

^KQQ 

O 

o 

^O 

O 

O 

O 

^DF 

o 

o 

^FF 

^KDF 

O 

^KDD 

o 

o 

^KDF 

^KDKL 

O 

O 

^KQQ 

o 

O 



O 

^KQKQ 

. K- 


R = diag 


[rd-  Rq*  Rq-  ^KD'  ^KQ  j 


and  jvij  = 0 for  all  Ly  except  aj2  " “''■21  * direct  phase 

variables  are  related  to  this  set  of  DQO  axis  variables  by  the  trans- 
formations 


-a be  ■ -DQO  '"a be 

where  a suitable  B~  ^ is  given  by 


B"  ^ V 


DQO 


1.9 


B" ^ = /2 /3 


eos  (i)t 


- sin  ujt 


1//2 


cos  (u)t  - 2tt/3)  - sin  ujt  - 2tt/3  1//’2 


cos  (u)t  + 2tt/3)  - sin  (ujt  + 2rr/3)  \IJ'2. 


1.  10 


Transforming  equation  1.8  to  direct  phase  quantities  via  1.9  gives 
® '"abc  = R B i^^jc  + A L B iabc  + L B i^bc  + L B — ^ 


or 


v^bc  = B-1  R B iabc  + [b’^  aL  B + B"^  L b]i 


*^abc 


1 

S '-®  — 


1.11 


1.  12 


Comparison  of  equations  1.12  and  1.5  indicate  that 
B-abc  = B"^  B^qo  B 

B-abc  = B"^  RdQO  ® 

As  all  possible  resistance  measurements  are  normally  made  in  direct 
phase  form,  no  transformation  is  necessary.  However,  inductance 
calculation  from  synchronous  and  transient  reactances  is  normally  in 
terms  of  Lqqq  and  matrix  entries  of  equation  1.4  are  related  to  Lqqq 


entries  as  follows: 
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1/6 


_Ldd  + LqQ  ] 


= ^KDKD 
= I"KQKQ 
= LfF 

= Maf  = LpD 

= ^KQQ 

= Lj^dF 

Of  course,  the  per  unit  denormalization  must  first  be  accounted 
for  to  acquire  direct  phase  parameters  which  represent  the  physical 

machine  and  reproduce  the  proper  level  voltage  and  current  variables. 

SECTION  IV 

LINEAR  MACHINE  SIMULATION 

The  parameter  set  in  fl]  was  used  to  simulate  a machine  in  the 
Sceptre  format.  The  parameters  were  unnormalized  to  give  the 
physical  phase  voltages  and  currents  under  normal  operating  conditions. 
The  model  and  parameters  of  Program  in  Appendix  A were  run  in 
various  load  and  switched  load  conditions.  Figures  1.3  and  1.4  show 
the  current  and  voltage  waveforms  on  the  A phase  as  the  open  circuit 
machine  has  a . 5jx  load  placed  on  the  A phase.  Figures  1.  5 and  1.  6 
show  the  result  on  the  A phase  variable  as  the  machine  has  a 3 phase 
line  to  neutral  . 5^  load  applied  on  phases  B and  C.  These 
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Lkkq 

Lff 

Mad 

Maq 
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computational  values  can  be  contrasted  with  the  experimental  results 
obtained  from  lab  tests  on  a similar  60  KVA  machine  at  WPAFB.  (Same 
model,  but  different  machine.)  These  results  are  shown  in  figure  1.7. 


The  open  circuit  machine  has  a 2 phase  line  to  neutral  load  of  approxi- 
mately . 5jt.  applied  on  phases  B and  C.  The  final  values  which  the 
machine  voltage  and  currents  approach  are  approximately  equal  for 
model  and  machine,  however,  the  time  constant  in  arriving  at  the  final 
values  is  slightly  different. 

This  indicates  that  field  coupling  and  resistive  values  are 
approximately  correct,  however,  damper  circuit  values  are  such  that 
their  action  is  too  dominant  on  the  machine  response. 

It  is  instructive  to  compare  the  simulation  time  for  this  machine 
in  the  Sceptre  format  with  results  utilizing  different  numerical 
techniques  for  the  time  solution  of  differential  equations.  For  problems 
involving  application  of  transient  load  networks,  the  most  severe 
accuracy  requirements  occur  immediately  following  the  transient,  thus 
figure  1.4  represents  a typical  loading  situation.  The  total  CPU  time 
on  an  IBM  370  Computer  for  problem  solution  excluding  Sceptre  set-up 
time  for  this  example  was  .22  minutes  or  13  seconds.  Typical  open 
circuit  to  short  circuit  computation  time  including  start  up  transient 
from  zero  initial  conditions  is  of  the  order  of  2 to  4 seconds /cycle 
for  this  linear  model. 
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Figure  1.7  60KVA  Machine  Voltage  Waveforms  - 25  Line  to 
Neutral  Short  (approx.  1 /2  J 


SECTION  V 

SIMULATION  OF  GENERATOR  SATURATION  CONDITIONS 


The  terms  in  the  machine  equation  1.  1 involving  the  time 
derivative  of  the  product  L^(t)  in  general  cannot  be  separated  into  the 
form  of  equation  1.5.  This  condition  arises  in  saturation  because  the 
inductance  elements  of  the  L matrix  are  fvinctions  of  both  time  and 
current.  However,  it  is  usual  practice  to  invoke  a separation 
principle  in  the  analysis  of  a saturated  synchronous  generator. 

This  separation  principle  [2]  is  based  upon  the  operation  of 

— r “1 

considering  the  derivative  I L(J_,  t)  Ht)  J as  the  sum 

d ' • ' 

L(i,  t)  TTT  i(t)  + L(i,  1)  i(t).  Under  this  assumption  it  can  be  seen  that 

5L  9i(t)  , j j rr.1 

a partial  derivative  of  the  form  Qi  ”§t  been  excluded.  The 

justification  of  this  saturation  principle  in  the  past  has  been  based  upon 
agreement  between  predicted  and  measured  results. 

The  approach  utilized  in  this  simulation  is  to  calculate  an 
equivalent  flux  with  a linear  combination  of  the  currents  of  the  windings. 
The  combination  of  currents  for  flux  calculations  is  taken  to  be 


^SAT  " 


whore 


Mgf  cos  Ojt 
Mgf  cos  (u)t  - 2tt/3) 
cos  (u)t  + 2tt/3) 


^FKD 
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It  is  seen  that  the  flux  calculation  is  such  that  the  linear  combination 
of  currents  produces  a rotating  magnetic  field  colinear  with  the  magnetic 
field  produced  by  the  field  winding.  The  inductance  values  are  modified 
as  a function  of  the  calculated  flux  and  a saturation  curve  for  the  machine. 

A saturation  curve  for  a 60  KVA  auxiliary  power  unit  is  shown  in 
figure  1.8.  This  curve  is  obtained  through  plotting  open  circuit  phase 
voltage  measurements  versus  field  current  when  rotating  at  synchronous 
speed.  The  zero  power  factor  curve  will  be  used  in  this  model  because 
of  the  absence  of  definable  power  factor  conditions  under  transient 
conditions.  If  the  saturation  curve  is  inverted  and  normalized  such 
than  an  equivalent  coupling  between  field  and  phase  is  computed  for  all 
values  along  the  saturation  curve,  the  result  shown  in  figure  1.9  is 
obtained.  Thus  the  inclusion  of  the  saturation  condition  is  equivalent 
to  modifying  the  machine  inductances  by  a ratio  determined  by  the 
calculated  flux  as  given  in  equation  1.  3 and  figure  1.  9. 

The  mechanics  of  including  this  effect  into  the  machine  model 
is  greatly  simplified  under  the  dependent  source  equation  format 
utilized  in  the  preceding  model.  A table  lookup  with  the  saturation 
curves  versus  calculated  flux  is  used  to  determine  the  reduction  factor 
for  the  inductance  values.  A circuit  containing  dummy  dependent 
sources  \^ich  calculate  the  flux  level  and  linkage  for  the  equivalent 
physical  orientation  of  the  phases  and  field  windings  is  introduced. 

The  source  value  for  the  physical  alignment  of  each  coil  with  the 
resulting  field  is  used  as  the  independent  variable  in  the  table  lookup. 


Neutral 
A.  C.  Voltage 


Figure  1,8  60KVA  Machine  Saturation  Curves 


Saturation 


I 


Figure  1.  9 Effective  Mutual  Inductance  (Field  to  Phase) 
Multiplier  for  Linear  and  Satviration  Regions 


The  resulting  multiplicative  factor  is  used  to  adjust  the  machine  inductive 
parameters  via  the  inductance  descriptive  equation.  The  program  for 

this  configuration  is  shown  in  Appendix  B. 

SECTION  VI 

SATliRATED  MODEL  RESULTS 

The  base  parameter  set  from  the  previous  linear  model  was 
used  in  this  simulation.  The  open  circuit  voltage  waveforms  and  phase 
voltages  and  currents  under  lOii  resistive  loads  produced  waveshaptes 
which  were  equivalent  to  the  linearized  model.  However,  the  results 
under  transient  loading  conditions  produced  waveshapes  with  high  3rd 
harmonic  content.  A typical  result  is  shown  in  figures  1.10  and  1.11.  In 

this  case  the  load  on  A phase  is  switched  from  10  a to  .01a.  The 
transitory  region  displays  the  same  general  damping  characteristic 
which  is  seen  in  the  linear  model.  The  steady  state  voltage  waveform 
under  high  current  (>100  amps/phase)  shows  the  characteristic 
rounding  effect  on  the  simusoidal  peaks  attributed  to  operation  in 
sa  turation. 

SECTION  VI 

CONCLUSIONS 

The  results  of  a specific  model  format  for  simulation  of  a 
synchronous  machine  have  been  presented  in  this  report.  This  forn^et 
satisfies  the  direct  phase  representation  necessary  for  the  inclusion  oi 
nonlinear  loading  networks  and  for  the  insertion  of  machine  nonlinear- 
ities. A method  of  insertion  of  the  saturation  effect  into  the  machine 
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Figure  1.11  Saturation  Model  - C Phase  Voltage  - Line  to 
Neutral  Short  (.01a)  on  A Phase 
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Figure  1.  10  Saturation  Model  - A Phase  Current  - Line  to 
Neutral  Short  (,01a)  on  A Phase 
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model  has  been  demonstrated  to  yield  realistic  results.  The  compu- 
tational time  of  both  linear  and  nonlinear  representation  is  such  that 
this  Sceptre  format  is  competitive  with  per/cycle  computation  times 
reported  in  the  literature.  The  waveform  results,  with  the  parameters 
used  in  this  simulation,  indicate  that  the  damping  characteristics  are 
extremely  short  lived  when  compared  to  the  experimental  results  of  a 
similar  a.p.u.  This  suggests  that  damping  circuit  equivalents  might 
be  better  described  through  simply  matching  the  waveform  damping 
characteristics  immediately  ‘‘ollowing  application  of  a load  transient. 
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SCEOTPF  SYSTEM  INPUT 


uHl  VEPSION 


S/3^'' 


OS  «FLSaS=  20.1 


CIRCUIT  CSSCRIPTION 
ELEMENTS 
Jl, 4,4-0=10.0 
R4L,''1-0  = 1P0'). 

RBL, 72-0=Tl 
RCL,73-0=T1 
R4,l-1 1=PRP 

RB, 2-12=P»P 

RC, 3-13=P»P 

033,13-23=  E0U4TI0N 
R22, 12-22=  ■^aU4fI0N 
Pll,  11-21=  ECU4TI0N 
RF,4-14=PRF 
ROtS-lSsPOfi 

Ra,6-16=PO0 


3  (-1.  ,PL42,-2.C<344» 
3 (-1.,PL42, 2.0944) 

3  (-1. ,PL42,0.P) 


E12,21-31=OQU4TICN  3 ( I L22  , PL42 ,-2  .0944 ) 
E13,31-41  = E0U.' TiON  3 ( I L 33  , P L4  2 , 2. 0944  | 
E14,41-51  = E0U1TICN  4 ( IL44,PM4F,0.C) 
E15,E1-61=  =QU4TinN  4 ( IL55  , PM4D  ,0 .0 ) 

E16, 61-71=  oyU4TlON  ( IL66,PM4Q,0.C» 

E21. 22-32=  ECt4TlON  3 ( I LI  1 , PL42 ,-2 .0944 ) 
E23, 32-42=  EQU4TIDN  3 (I L 33 , PL 42 , 0. 0 ) 


E24,42-F2= 

625.52- 62= 

026.62- 72= 

631.23- 33= 

032.33- 43= 

634,43-53= 

635.53- 63= 

636.63- 73= 

641 . 14- 24= 

601.15- 20= 

642.24- 34= 

652.25- 35= 

053,30-0  = 

643.34- 44= 

661.16- 26= 

662.26- 36= 


62U4T10N 

6QU4TinN 

6CU4TI0N 

EQU4T10N 

OQU4TICN 

6QU4T10N 

FaU4TI0N 

EQU4TinN 
EQUATION 
OCUATION 
EQUATION 
0QU4T  ION 
EQUATION 
oQUATION 
OQLATION 
EQUATION 


063,36-0  = EQUATION 
L11,C-1=  OQUATION  1 
L22,0-2=  OQU4TI0N  1 
133,9-3=  EQUATION  1 


4 ( IL44,PMAF,-2.0944) 

4  ( IL55,PMAD,-2.0944) 
o ( IL66,PWAQ,-2.0944) 
3 HLl  1,PLA2  ,2.0944) 

3 ( IL22,PLA2,9.0) 

4 ( IL44,PMAF,2.0944) 

4 ( IL55,PMAD,2.0944) 

5 ( IL66,P»'A0,2.0944) 

4 ( ILl 1,PYAF ,0.0) 

4 < ILll,P5'A0,0.O) 

4 ( IL22,P«AF ,-2.0944) 
4 (IL22,PMAD,-2.0944) 
4 ( 1L33,PMAC,2.0944) 

4 (IL33,PMAF,2.0944) 

5 (ILll,PMAC,P.O) 

5  <1L22,P»'A0,-2.0944) 
5 ( IL33,PMAO,2.0944) 
(PLA0,PLA2,0.0) 
(PLA0,PLA2,2.O944) 

< PLAO, PLA2  ,-2.0944  ) 


L44,r-4=PLO0 

L55',C-5  = PLKK0 
L66,9-6=PLKKQ 
f»45  ,L44-L55  = PPFKC 


►«12,L11-L22  = 
•'13,L11-L33= 
P15,L1 1-L55= 

^14,L1 1-L44= 
)»23,L22-L33  = 

M24,L22-L44= 
H25,L22-L50= 
M34,L33-L44= 
P35,L33-L55= 
P16,L1 1-166= 
)*26,L22-L66  = 
P36.L33-L66= 


.OQUATION  1 
EQUATION  1 
EQUATION  9 
EQUATION  9 
FQUATION  1 

equation  9 

EQUATION  9 
OQUATION  9 
EQUATION  9 
FQLATION  B 
EQUATION  B 
OQUATION  B 


(PMS0,PLA2,-2.0944) 
(P*^SC,PL  A2,2.0944) 

( 1.0,PMA0,0.O) 

( l.C,PMAF ,0.0) 
(P«S0,PLA2,0.O) 

I 1.0,PNAF,-2.0944) 
( 1.0,PMAr),-2.0944) 

( 1.0,PMAF,2.0944) 

( 1.0,PMAD,2.C944) 

I-l.O, PMAQ.0.0  ) 
(-1.0,PMAQ,-2.0944) 
(-1 .0,PPAO,»2.0944  ) 
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DEFIN^r  PfioauFTFRS 
P®C  = 

PPO= 

P^P  = .''<>57 
PRF  = .''9‘'' 

PL  A'>  = . 

®L  A ? =.  r'''“rjo 
PLF‘^  = .CS'^ 

PLKKO-  “T*-» 

PLKKC=  ."'■'T51'5 
P'*S''  = -.0C'>2?1 
p<*FKn=  .or^'fl 

PM/if  =-  ,nr<)0  a 
PM4C=-.CC04 
o*^AC=-  .0014  ) 

PW=?C12. 

FIJNCtioms 

FauATTc^'  1 (A,e,c»  = 
FJU1TTC»'J  ^ «A,«,C)  = 
cqUATton  4 ( A, B,CI  = 
FQLiATI^^,  15  (A,H,C»  = 
EQUATTOM  (A,P,C)  = 
FjuATirM  B (A,P,CI  = 
FQUAticn  9 (A,P,C(  = 

T1 

r.r.  IF'':''.  ,r).9oc  ,xo-3c. 


( At3*0C0S(2.4PW*TIMF  + C> I 
(2.*A*P*Pw*DSIN(2.*PW'‘TIME-i-C)  » 
(A*P«PW*r>SIN(PW*TIMF  + CM 
( A*9* PW*DCnS(PW*TI“E4C) J 
( 2.*A*8*0SIN( 2.*PW*TIME*C) ) 

( A*P*DSIN<PW*TIME*Cn 
( A»B*DCnS(  PW*TI*1E  + Cn 

9091, .5, .02, .5 


CUTPUTS 

V»AL ,VPPL , ILll , ILf6,PLCT 
O'IN  CONTROLS 

intrcoa-ficn  p-^ut ine  = impl ic it 

STOP  ’'I«F  =0,02 


CNC 
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APPENDIX  B 


SCEPTRE  SATURATION  MOUEL  PROGRAM 


II 


SCEPTRE  SYSTEM  INPUT 


4/72  VERSION  — S/T#.0  — nS  RFLFASE  20.1 


CIRCUIT  DESCRIPTION 
ELEMENTS 

ES1.0-A^SATnLllfIL22t  U33.IL44,  IL55, 0.0, TIME  » 
ES2,B-C=ESATI1L11 ,1L22 , 1 L33 , 1 L44, 1LS5, -2.0944 ,T1 MF » 
ES3,A-B=FSAT<1L11  ,IL22,IL33,1L44,1L5S,^2.0O44,’^Tmfi 
ES4,C-0=PSAT(IL11,IL22,1L33,IL44,IL55,1.0,TTMF) 
RDEL,0-Ci=lQOO. 

Jl, 44 -0=10.0 
RAL,71-0=T2 
RBL, 72-0=10. 


RCL,73-0=10. 

RA, 1-11=PRP 

RB, 2-12=PRP 

RC, 3-13=PRP 

Rll, 11-21*  EQUATION 
R22, 12-22*  EQUATION 
R33, 13-23*  EQUATION 
RF,4-14=PRF 
R0,5-15=PR0 
R0,6-16=PRQ 


3 «-l.,PLA2, 0.0, TABLE  1 (FSIM 
3 1-1.  ,PLA2, 2. 0944, TABLE  1 ♦FS2H 
3 (-1.  ,PLA2, -2. 0944, TABLE  1 »FS3U 


E12 

E13 

E14 

E15 

E16 

E21 

E23 

E24 

E25 

E31 

E26 

E32 

E34 

E35 

E36 

E41 

E42 

'43 

..63 

E62 

E61 

E53 

E52 

ESI 

L44 

L55 

L66 

M45 

L33 

L22 

Lll 

M16 

M15 

M14 

M13 

M23 


3 UL22,PLA2, -2. 0944, TABLE  1 lESll) 

3 (IL33,PLA2, 2. 0944, TABLE  1 fESTM 

4 nL44,PMAE, 0.0, TABLE  1 IFSIM 
I IL55,PMAD,0.0,TABLE  1 (ESm 
«1L66,PMA0, 0.0, TABLE  1 (ESIH 
«1L11,PLA2, -2. 0944, TABLE  I IPS?)) 
ML33,PLA2, 0.0, TABLE  1 »FS2H 

42-52=  EQUATION  4 «1 L44,PMAE , -2 .0944, TABL E 1 (ES2M 
52-62=  EQUATION  4 M L55,PMAD, -2. 0944, TABLE  1 <FS2)) 
«1L11,PLA2, 2.0944, TABLE  1 (ES3U 
UL66,PMA0, -2. 0944, TABLE  1 IES7»> 
<1L22,PLA2,0.0,TABLE  1 «ES3n 
UL44,PMAE, 2.0944, TABLE  1 »ES3|» 
nL55,PMA0,2.0944,TABLF  1 ftSBU 
( 1L66,PMA0, 2. 0944, table  ! (PSB»| 
(ILll.PMAE, 0.0, table  1 (ES4»» 
UL22,PMAE, -2. 0944, TABLE  1 IES4)1 
UL33,PMAE, 2. 0944, TABLE  1 (ES41) 

( IL33,PMA0, 2.0944, TABLE  1 (ES3») 
aL22,PMA0,-2.0944,TABLF  1 »ES?n 
<IL11,PMAQ, 0.0, TABLE  1 |FSn» 
<1L33,PMAD, 2. 0944, TABLE  1 IES4J) 
•IL22,PMAD, -2. 0944, TABLE  1 |FS4)1 
(IL11,PMA0,0.0,TABLF  1 (FS4I) 

0-4=  EQUATION  6 IPLEE, TABLE  1 1ES4M 
0-5=  EQUATION  6 IPLKKD, TABLE  1 <ES4U 
0-6  =PLKKQ 

L44-L55*  EQUATION  6 IPMEBO, TABLE  1 «FS4n 
0-3=  EQUATION  1 I PL AO, PL A2, -2 .0944 , TABLE  1 IES3)) 
0-2=  EQUATION  1 |PL AO, PL A 2 , 2.0944 , TABLE  1 (ES2JI 
0-1=  EQUATION  1 IPLA0,PLA2, 0.0, TABLE  1 lESTU 
L1I-L66=  EQUATION  S «-l  .0,PMAQ  ,0.0, TABLE  1 FESIU 

« i.o,pmad,o.o,tabip  1 |F?4n 

« I.0,PMAE,0.0,TABLE  1 |FS4»» 
»PMS0,PLA2, 2. 0944, TABLE  1 (ESI)) 
(PMS0.PLA2, 0.0, TABLE  1 (FS2)) 


21- 31=EOUAT10N 

31- 41 =E0UAT10N 
41-51=EQUA710N 
51-61*  EQUATION 
61-71=  EQUATION 

22- 32=  EQUATION 

32- 42*  EQUATlLiN 


23- 33=  EQUATION 

62- 72*  EQUATION 

33- 43=  EQUATION 
43-53=  EQUATION 
53-63*  EQUATION 

63- 73=  EQUATION 

14- 24=  EQUATION 

24- 34*  EQUATION 

34- 44=  EQUATION 
36-0  e EQUATION 
26-36=  EQUATION 
16-26=  EQUATION 

35- 0  = EQUATION 

25- 35*  EQUATION 

15- 25=  EQUATION 


L11-L55*  EQUATION  9 
L11-L44*  EQUATION  9 
L11-L33*  EQUATION  1 
L22-L33*  EQUATION  I 
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H2i,L22-t55=t  EOUATION  9 
H26,L22-t66*  FeUATION  8 
M35,L33-t55*  EOUATION  9 
M36tL33-iA<>*  EQUATION  8 
M12,un-L22*  EQUATION  1 
N24,t22-i.A4=  EQUATION  9 
N34,L33-CA4*  EQUATION  9 
OEFINEO  PARAMETERS 
PR0».077? 

PR D=. 0778 
PRP*.04f,7 
PR F*. 0905 
PLAO* .00044? 

PLA2*.0CO09 
PLFF».085 
PLRKO*. 00078 
PLKKO*. 000515 
PHSO*-. 000221 
PMFK0  = .«)78 
PMAF*-. 00608 
PMAO*-.0004 


« I.O,PMAD,-2.0944,TAPlF  1 
|-l.<»tPMA0.-2.0944,TAPLP  ! 
I 1.0, PNAO,?. 0944, TABIF  1 

1-1  .0,PM*0,42.0944,7APLP  1 
(PNS0,PLA2,-2.0944,TAALF  ? 
( 1.0,PMAF,-2.0944,TA«LF  1 
< 1.0,PNAF,2.0944,TA»L»-  1 


( FS4  n 
»FS?»» 
»FS4)I 
|PS3n 
ffsin 
(FS4)» 
EPS4M 


PH AD»-. 00060 
PW*2512. 

FUNCTIONS 

EQUATION  9 (A,B,C,0»  = H A*0*DCOS <PW*TINE4C II*" I 
EQUATION  8 (A,8,C,0)  = ((  A*BAOSIN ( PW*T1 HE*C » »*0) 
EQUATION  7 (A,B,C,D»  * ((  2.PAPB40SINI 2.*PHPT1mF4C  11*01 
EQUATION  6 IA,B)  = (A*BI 

EQUATION  5 (A,P,C,0)  * MA*n*PW*OCCSIPM*TlNF*Cl )*0» 
EQUATION  4 <A,B,C,0)  = ( (A*B*PW*0S1N  I PW*T  I“E*r  ))••>) 
EQUATION  3 (A,B,C,0l  = n2.*A*B*PW*0SIN1 2.*“M*TTHr»C) )*0) 
EQUATION  I (A,b,C,D»  = H A*B*DCOS( 2. *PW*T IMP*C 1 1*0) 

TABLE  I 


0.0,1  .0 ,1 .0,  1.0, 1.5, .8,2 .0 , .7, 2. 5, .54, 3.0, .44 ,4.0,  . 36 , 5 .0  , . 3 

T2 

0.0,10.0,0.009,10.0,0.0091 ,0.01 ,0.02,0.01 
OUTPUTS 

1L11,VRBL,VRCL,IL55,IL66,PLUT 
RUN  CONTROLS 

INTEGRATION  R0UT?NE*1NPLIC IT 
STOP  TINE  =0.03 
MINIMUM  STEP  SIZE  * l.OE-11 
END 

INDICES  FOR  state  VARIABLE  DERIVATIVE  NAMES 


1 OL55 

2 DL66 

3 DL33 

4 OL22 

5 DLIl 
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IV  G L*=VEL  20.1 


FSAT 


DATE  s 76195 


12/Al/<>6 


DUUBLF  PRECISION  FUNCTION  FS AT < AL t BL, CL,OL t E L »FL » T INE ) 

IMPLICIT  REAL*8( A-J,L-M,o-7) , INTCGFR*4 IK ,N ) 

PW=2512. 

PMAF=-.00608 

PLFF=.085 

PVFKr>=.0078 

FLUX=PMAF*AL*OCUS ( PH*TI MF  > ♦pmaf*BL*OCOS ( PW^TIME-a .0944) ♦PNAF*Ct* 
* nCCSIPW*TIME-*’2.0944)-»PLFF*nU  + «>MFKD*EL 
FSAT=r)ABSIUCOS(PW*TIMF-FL  )*FLUX» 

IF  IFL.FO.l.O)  FSAT=OARS(FLUX) 

RETURN 

END 


! 
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